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ABSTRACT 

We present a detailed study of the early phases of the peculiar supernova 2011ay 
based on BVRI photometry obtained at Konkoly Observatory, Hungary, and optical 
spectra taken with the Hobby-Eberly Telescope at McDonald Observatory, Texas. The 
spectral analysis carried out with SYN++ and SYNAPPS confirms that SN 2011ay 
belongs to the recently defined class of SNe lax, which is also supported by the prop¬ 
erties of its light and color curves. The estimated photospheric temperature around 
maximum light, Tphot ~ 8,000 K, is lower than in most Type la SNe, which results in 
the appearance of strong Fell features in the spectra of SN 2011ay, even during the 
early phases. We also show that strong blending with metal features (those of Till, 
Fell, Coll) makes the direct analysis of the broad spectral features very difficult, and 
this may be true for all SNe lax. We find two alternative spectrum models that both 
describe the observed spectra adequately, but their photospheric velocities differ by at 
least ~ 3,000 km s“^. The quasi-bolometric light curve of SN 2011ay has been assem¬ 
bled by integrating the UV-optical spectral energy distributions. Fitting a modified 
Arnett-model to Lyoiit), the moment of explosion and other physical parameters, i.e. 
the rise time to maximum, the ®®Ni mass and the total ejecta mass are estimated as 
trise ~14±1 days, ~0.22±0.01 Mq and Mgj ~ 0.8 Mq, respectively. 

Key words: supernovae: general - supernovae: individual (SN 201 lay) 


1 INTRODUCTION 

Type lax supernovae (SNe lax), previously labeled as 
SN 2002cx-like ones, are spectroscopically similar to SNe 
la, but thought to have lo wer velocities at m aximum light 
and lower peak magnitudes dFolev et al.ll20lJ l. There are 25 
members of this class known to date, but it is suggested that 
in a given volume there may be 30-40 SNe lax for every 100 
SNe la. 

The major question about SNe lax is the nature of their 
progenitors. The similarity of their light curves to those of 
“normal” SNe la (especially their rise times), their spectral 
features, and the lack of X-ray and radio detections sug¬ 
gest that they originate from compact progenitors, probably 
white dwarfs (WDs). At the same time, explosions causing 


the total disruption of WDs (either via detonation or via 
pure deflagration) are not in agreement with the observa- 
tional results. T his sta tement was recently strengthened by 
IMcCuUv et ahl ll2014ah who analyzed late-time spectra of 
two SNe lax. 


Partial disruption of WDs seems to be a feasible 
model to expl ain the properties of this new class o f stel¬ 
lar explosions dFolev et al.lf20lJ : lMcCullv et al.ll2014al l. This 
model seems to give a plausible explanation for the di¬ 
versity of this class (including extremely faint and low- 
velocity explosions, like SN 2008ha). It is also supported 
by recently published results of numeric al 3D simulations 


( Jordan et al.1 20121: Kromer et al 


[ 201 ^. As iFolev et al.l 
1 2 OI 3 II and Wang. Justham fc Han (1201311 proposed, WDs 
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accreting material from their helium-star companion in bi¬ 
nary systems may be good candidates for the progenitors 
of SNe lax. This assumption has been recently strength¬ 
ened by the detection of the lumi nuous, blue progenito r 
system of the type lax SN 2012Z llMcCullv et al.l f2014bll . 
At the same time, the observed properties of this object, 
and maybe those of other SNe lax, can be also explained 
with the pulsational d elayed detonation of a single WD (see 
IStritzinger et al]|2015l . and references therein). 

Other ideas, e.g. the theory of the “.la” events (which 
are thought to be luminous helium shell flashes of WD s in 
binary systems with low accretion rates. iKldstenllioOTl l. do 
not se em to give self-consiste nt description of the observa¬ 
tions dFolev et al.ll200^ . l2013l l. Note, however, that the .la 
models might be good alternatives to explain the extremely 
low-luminosity events, like SN 2008ha. 

On the other hand, so far almost all SNe lax have 
been found in late- type galaxies with recent, strong sta r 
formation activity jLvmau et al.l l2013l : iFolev et al.l l2013ll . 
This fact might lead to the suspicion that SNe lax actually 
could be low-luminosity core-collapse ( CC) events emerging 
from 7-9 Mp), stri pped-envelope stars (IValenti et al.l [ 2 OO 9 I : 
iLvman et al.ll201^ 1. even though the observed properties of 
SNe lax are more similar to those of la than CC SNe. In this 
scenario the explosion might not have enough energy to eject 
the entire stellar envelope , thus, part of the ejecta falls back 
to th e central remnant (IValenti et al.l 120091 : iMoriva et al.l 
|201(J) . Note, however, that SN 2008ge, another SN lax, is 
situated in a region that has no detectable star-formation, 
thus, there are probably no massive stars near the SN site 
llFolev et al.ll2010bl. l2013l'l . 


To understand the general properties of SNe lax, 
it is important to study single, well-observed ob¬ 
jects as thoroughly as possible. Up to now, de¬ 
tailed photometric and spectroscopic analyses have 
been published only on a few of them: SN 2002cx 


llLi et al.l l2003l: Branch et al. 

20041: Jha et al. l2006l). 

SN 2005hk llChornock et al. 

20061: IphilliDS et all l2007l: 

ISahu et al. 2008: 

McCullv et al.l l2014ah. SN 2007ad 

dMcClelland et al. 2010f), SN 2008A (iMcCullv et al.l 2014311. 

SN 2008Ee llFolev et al. l2010bl). SN 2008ha IFolev et al.l 

I 2 OO 9 I: Valenti et al. 

20091: Folev et all 2010a), SN 2009ku 

iNaravan et al. 2011 

), SN 2010ae IStritzinger et al. l2014h. 

SN 2012Z (IMcCuIIv et al.l l2014bl: IStritzinger et al.l 20151: 


Moreover, Folev et al.l (l2013l j presented light curves and 
spectra of some other SNe lax for which detailed studies 
are in progress. 

Here we present a detailed analysis of SN 201 lay, an¬ 
other SN lax. This object was discovered on March 18.18 UT 
by the KAIT/LOSS pr ogram with an appare nt, unfiltered 
brightness of 17.7 mag (iBlanchard et alTl201ll '). The SN lo¬ 
cated 9.3”east and 1.4"south from the center of NGC 2315. 
The average Hubble-flow distance of the host galaxy is 
86.9±6.9 Mpc, as listed in the NASA/IPAC Extragalac- 
tic Database ( NED). Based on a quick analysis of an 
early spectrum, IPogge. Garnavich fc Pedanil ll201lll classi - 
fied SN 2011ay as a peculiar la. while lSilverman et al.l (1201 il l 
reported that it m ay belong to the SN 2002cx-like events. 
iFolev et al.l ll2013h confirmed the classification and also pub¬ 
lished some optical photometry and spectra of this object. 

Our motivation was to carry out a detailed analysis us- 
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Figure 1. BVI composite image of the SN 201 lay and its host 
galaxy on 2011 March 24 (0.6/0.9 m Schmidt-telescope, Konkoly 
Observatory, Hungary). The numbered objects were used as local 
comparison stars. 


ing data obtained during the early photospheric phase of the 
SN. We present our simultaneous photometric and spectro¬ 
scopic observations in Section [S] while the light and color 
curves, the UV-optical spectral energy distributions (SEDs) 
and the results of modeling the spectra and the light curve 
are shown in Section O In Section |T] we discuss our results 
and present our conclusions. 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Photometry 

Ground-based photometric observations for SN 2011ay were 
obtained from the Piszkesteto Mountain Station of Konkoly 
Observatory, Hungary. We used the 0.6/0.9 m Schmidt- 
telescope and the 1.0m RCC-telescope, both equipped with 
Bessell BVRI filters. Figure [T] shows a BVI composite image 
of the SN and its host galaxy. 

We used the daophot/allstar task in IRAfQ to carry 
out PSF-photometry on the SN and two nearby compari¬ 
son stars (see Figure [TJ. Because the SN appeared close to 
the center of the nearly edge-on host galaxy, the value of 
background flux level was estimated interactively in every 
case. The background level estimates produced by the fit- 
skypars task were adjusted by eye and the background was 
subtracted manually for each frame. 

The instrumental magnitudes were transformed to the 
standard system applying the following equations: 


V - n = Cv • (U - /) + Cu 
{B-V) = Cbv ■{b-v) + Cbv 
{V — R) = CvR • (u — r) -b (^vR 

iV -I) = Cvi-{v-i) + Cvi, (1) 

where lowercase and uppercase letters denote instrumen¬ 
tal and standard magnitudes, respectively. The color terms 
(Cx), determined by measuring Landolt standard stars 
during photometric conditions close to the epochs of the 


^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. BVRI magnitudes of SN 201 lay from the Konkoly Ob¬ 
servatory, Hungary 


JD“ 

B (mag) 

V (mag) 

R (mag) 

I (mag) 

5645.3 

16.89(.07) 

16.97(.05) 

16.75(.06) 

16.74(.10) 

5647.3 

16.85002) 

16.71(.05) 

16.52006) 

16.33(.08) 

5649.3 

16.99003) 

16.62(.07) 

16.44006) 

16.31008) 

5650.3 

16.99(.06) 

16.54004) 

16.50005) 

16.28(.08) 

5651.3 

16.97005) 

16.65007) 

16.42006) 

16.21(.07) 

5655.3 

17.30(.06) 

16.62007) 

16.31006) 

16.10007) 

5657.3 

17.38(.07) 

16.63(.06) 

16.28006) 

16.09(.07) 

5661.6 

18.02(.15) 

16.86006) 

16.36005) 

16.21(.07) 

5669.3 

18.82(.23) 

17.71003) 

16.88004) 

16.43002) 

5671.3 

19.05(.12) 

17.78009) 

17.05(.03) 

16.48(.02) 


Notes. JD—2,450,000. Errors are given in parentheses. 
Table 2. Swift magnitudes of SN 2011ay 


JD“ 

UVWl (mag) 

U (mag) 

B (mag) 

5645.5 

17.70(.12) 

16.53(.07) 

16.59(.05) 

5648.0 

17.98011) 

16.63(.07) 

16.66(.05) 

5650.1 

18.22013) 

16.75(.06) 

16.67005) 

5651.8 

18.32014) 

16.91(.08) 

16.80005) 

5653.5 

18.17013) 

17.14(.10) 

16.79(.06) 

5655.6 

18.66025) 

17.26(.12) 

17.01(.08) 

5657.6 

18.73019) 

17.36(.10) 

17.18007) 

5660.1 

18.83022) 

17.88018) 

17.54011) 


Notes. 2,450,000. Errors are given in parentheses. 

SN observations, were the followings: Cv = 0.047, Cbv 
= 1.228, CvR ~ 0.960, Cvi = 0.934. Because simultane¬ 
ous observation of the Landolt field and the SN field was 
not possible, we applied Eq. o to the differential magni¬ 
tudes Am = msN — mci ,2 constructed from the instru¬ 
mental magnitudes of the SN and the comparison stars, 
thus, eliminating the zero points (Cx)- Then, using the 
standard BVRI magnitudes of the comparison stars calcu¬ 
lated f rom their SDSS-magnitudes ( based on the calibra¬ 
tion of ijordi. Grebel fc AmmMl2006l l. the standard magni¬ 
tudes of SN 2011ay were obtained in the Johnson-Cousins 
system. The results are presented in Table [T] Errors (given 
in parentheses) contain both the uncertainties of the PSF- 
photometry and the standard transformation. 

The optical data were supplemented by the available 
Swift/UVOT data (reduced using standard HEAsoft tasks). 
Individual frames were summed with the uvotimsum task, 
and magnitudes were determined via aperture photometry 
using the task uvotsource. The object could be roughly sep¬ 
arated from the galactic nucleus only in channels b, u, and 
uwl, while it got too faint after ~2 weeks. The results of 
Swift-photometry are shown in Table (2] 

2.2 Spectroscopy 

Optical spectra of SN 2011ay were obtained with the 9.2m 
Hobby-Eberly Teles cope (HET) Ma rcario Low Resolution 
Spectrograph (LRS, iHill et ahl Il998l l at McDonald Obser¬ 
vatory, Texas, between March 24 and April 19, 2011. LRS 


Table 3. Log of spectral observations obtained with HET LRS 


Date 

JD“ 

Phase^ (d) 

Exp. time (s) 

2011-03-23 

5643.2 

-10 

1500 

2011-03-24 

5644.2 

-9 

1800 

2011-03-26 

5646.2 

-7 

1500 

2011-03-27 

5647.2 

-6 

1500 

2011-03-28 

5648.2 

-5 

1500 

2011-03-29 

5649.2 

-4 

1500 

2011-03-31 

5651.2 

-2 

1500 

2011-04-14 

5665.1 

+12 

1500 

2011-04-19 

5670.1 

+n 

1500 


Notes. OljD—2,450,000; with respect to the moment of max¬ 
imum light in V-band. 



Figure 2. HET spectra of SN 2011ay. Epochs are given relative 
to V-band maximum. The spectra are shifted vertically for better 
visibility. 

has a spectral coverage of 4,200-10,200 A and a resolving 
power of A/AA 600. The data were reduced with standard 
IRAF routines. Table |3] contains the journal of the spectro¬ 
scopic observations, while the extracted, wavelength- and 
flux-calibrated spectra are collected in Figure [2] 


3 DATA ANALYSIS AND RESULTS 
3.1 Light curves and colors 

The light curves consisting of standard BVRI and Swift 
magnitudes are shown in Figure [S] The Konkoly B and V 
curves are general ly consistent with the ones published by 
iFolev et al.l ll2013l L but the sampling of our data around 
maximum is better (at longer wavelengths they observed 
SN 20011ay with r an d i filters). As seen in Fig. |3] the data 
by iFolev et al.l ll2013ll do not provide new information in 
addition to our photometry during the early phases. While 
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Figure 3. Standard BVRI (Konkoly Observatory) and Swift light 
curves of SN 2011ay (not corrected for reddening). Th e open cir¬ 
cles denote the B and V data from l^lev et aO J2013h . 
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Figure 4. Color curves of SN 2011ay (corrected for Milky Way 
reddening, see the text for details). The values plotted with 
open circles were calculated using the B and V datapoints of 
iFolev et al.l |2oT3). 


they followed the SN longer, their late-time photometry is 
much more uncertain, probably due to the increasing con¬ 
tribution from the host galaxy background flux. Thus, for 
further analysis we used only our data, but n ote that the 
result s are consistent with the photometry by iFolev et al.l 
11201311 . 

By fitting low-order polynomials to the data around 
the maxima, the light-curve parameters could be deter¬ 
mined with relatively low uncertainties. These parameters 
are listed in Table 01 Note that for normal SNe la the mo¬ 
ment of B maximum is usually used as a reference epoch. 
Instead, we used the time of V maximum, which is closer to 
the moment of maximum luminosity, to be consistent with 
the analysis presented by IFolev et al.l (l2013ll . 

The color evolution of SN 2011ay (Figure|4]) is very sim¬ 


ilar to other SNe lax (see e.g. El et ahllioosl : IPhillips et al.l 
l2007l : lrblev et al.ll2013fl : it gets redder after V -maximum un¬ 
til it reaches a nearly constant color at 15-20 days after max¬ 
imum. Although the interstellar matter (ISM) in the host 
galaxy may contribute to the total reddening, it is proba¬ 
bly low, because the SN appeared very blue, {B — V) <0 
mag, at the earliest observed epochs. Thus, the color curves 
have been corrected only for Milky Way reddening, using 
E{B — V) — 0.081 r nag. Ay = 0.23 mag, Ar = 0.18 mag, 
and Ai = 0.13 mag (ISchlaflv fc Finkbeinerllioilll . 


3.2 Spectroscopic analysis 

In the following, we present detailed analysis of some of the 
HET spectra shown in Fig. [2] and Table [S] The —9d, —2d 
and +17d spectr^ were selected for this purpose, because 
they have the highest signal-to-noise and show the largest 
differences in their spectral features, i.e. the most notice- 
abl e evolution. Moreo ver, two additional spectra obtained 
by IFolev et ahl ll2013ll were also included in the sample of 
analyzed spectra. These two spectra were taken with the 
Lick/Kast spectrograph on April 2 and 5, 2011, at Od and 
+3d phases, respectively. 

T he application of the Su pernova Identification (SNID) 
code dBlondin &: Tonr\1 l2007l ~l to the pre-maximum spec¬ 
tra of SN 2011ay revealed some similarity with spectra of 
normal Ic SNe. However, using SNID for the +17d HET 
spectrum we found that the relatively weak Si, dominating 
Fell and strong Call features clearly imp ly that SN 2011a y 
most closely resembles SN 2005hk (see e.g. lSahu et al.lf2008l ~l. 
i.e. it belongs t o the lax s ubclass, as f i rst po inted out by 
ISilverman et al.l (l201l[l and IFolev et al.1 (l2013fl . 

For more quantitative an alysis, we applied the 
SYN++ and SYNAPPS codes (iThomas. Nugent fc Me^ 
I2OIII) , which are based on the original SYNOW code 
i Jefferv fc Branch! Il990l : iHatano et al.l[l999l l. to model the 
selected spectra in order to reveal the chemical composition 
and other physical properties of the ejecta of SN 2011ay. 
SYNOW/SYN++/SYNAPPS assumes a fully opaque pho¬ 
tosphere located at a specific velocity v^oh ot, where the 
electr on-scattering optical depth Tea = 1 (iHatano et al.l 
Il999ll . in homologously expanding ejecta. The spectral fea¬ 
tures are assumed to be due to pure resonant scattering 
of the photons originating from the (assumed) blackbody 
emission of the photosphere. The scattering of photons take 
place in the partly transparent ejecta above the photosphere, 
where the strong velocity gradient (due to the rapid, homolo¬ 
gous expansion) makes the scattering regions appear as thin 
planes perpendicular to the line of sight for each photon 
frequency/wavelength. This Sobolev-approximation, among 
others, is clearly a limitation for the applicability of such a 
code. Nevertheless, it was found to be surprisingly effective 
and useful for modeling P Cygni line prohles and identify¬ 
ing SN features (se e e.g. IParrent, Friesen fc Parthasaratlm 
l20l4IParrentll2014l'l . 

Before the fitting, all spectra were corrected to Milky- 
Way extinction (adopting E{B — V) = 0.081 mag as 


^ We use the moment of V-band maximum light for assigning 
phases to the observed spectra, in accord with IFolev et all ll2013l i 
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Table 4. Light-curve parameters of SN 2011ay 



B 

V 

R 

I 

^max 

(JD-2,450,000) 

5646.6(1.1) 

5653.6(0.4) 

5657.2(1.1) 

5656.2(0.8) 

Peak mag. 

16.89(.08) 

16.56(.08) 

16.30( .07) 

16.08(.08) 

Peak abs. mag. 

-18.15(.17) 

-18.39(.18) 

-18.60(.17) 

-18.76(.18) 

Amis (mag) 

1.11(.16) 

0.95(.08) 

0.82(.07) 

0.40(.08) 


Notes. Errors are given in parentheses. 


above) and the redsh ift of the host galaxy using 2: = 0.021 
dMiller fc Owenll200ll) . 

Since parametrized models, like those from 
SYNAPPS/SYN++, depend heavily on the a-priori 
identihcation of ions and key parameters such as the 
velocity at the photosphere, we built two alternative models 
based on slightly different assumptions. 

In the first model (hereafter Model A) all featnres were 
assumed as photospheric, i.e. the line formation region for 
all ions extends from the top of the ejecta, parametrized as a 
fixed Vmax = 40, 000 km s”*^ at all phases, down to the sharp 
edge of the photosphere, given as Vphot in velocity coordi¬ 
nates. The density structnre of the ions distributed within 
the line forming region is treated as a simple exponential, 
defined directly for the optical depth of a given feature as 

t{v) = T(Ure/) X exp[-(u - Ure/)/ne], (2) 

where v^ef is a fixed reference velocity (arbitrary, but set 
close to Vphot) and is the e-folding width of the optical 
depth proHle. SYN-I—I- uses a quasi-LTE approximation for 
computing the relative strengths of the features for the same 
ion: the model specifies a single T(vref) optical depth for a 
pre -selected “reference ” feature of the given atom/ion (see 
e.g. lHatano et af]ll999l . for the list of reference lines for every 
ion), then the optical depths of all other same-ion features 
are computed assuming Boltzmann-excitation. The excita¬ 
tion temperature T^xc can be set differently for every ion, 
mimicking a non-LTE-like excitation, although it is clearly 
very far from a true, self-consistent non-LTE treatment of 
the problem. Thus, in Model A, all features are assumed to 
be formed down to the photosphere, each ion can have a 
different excitation temperature (Toxc), a different reference 
line optical depth (r(Ure/)) and different scale height of its 
line forming region (ve). 

The alternative model (hereafter Model B) was built 
by assuming two essential differences from Model A: in¬ 
stead of requiring that all features start to form at Vphot, 
the minimum velocity of their line-forming region, Vmin, 
may be at higher velocities than Vphot- S uch features, hav¬ 
ing v -min > Vphot, are called “detached” jjefferv &: Branchl 
Il99d l. This model may be more flexible than Model A, be¬ 
cause Vmin can be different for every ion, but also less con¬ 
strained, because the increasing number of parameters may 
decrease the uniqueness of the best-£t model. In addition, 
since Vphot may become much less constrained in this model, 
we set th e initial value of Vphot at ~ 6,000 km s~^ as esti- 
mated bv ISilverman et ^ |20l3) an d iFolev et al.l (l2013l j for 
SN 2011ay particularly from the absorption minimum of the 
feature around ~ 6200 A, which was thought to be due to 
Sill (but see Sect l3.3l for a discussion of this issue). 


Table 5. The photospheric velocities and temperatures for Model 
A and B, as found by SYNAPPS 


Date 

Epoch 

(days) 

phot 

(km s“l) 

/-pA 

phot 

(K) 

phot 

(km s~^) 

p'B 

phot 

(K) 

March 24 

-9 

9830 

8710 

6500 

10,000 

March 31 

-2 

9640 

7460 

6500 

10,000 

April 02 

0 

9280 

7260 

6000 

8760 

April 05 

-b3 

9090 

6000 

6500 

8620 

April 19 

-bl7 

8870 

5300 

4100 

6490 


The chemical composition of both models was assem¬ 
bled based on the major feat ures found in the s pectra of 


other SNe lax, namely 2002cx llBranch et al. 

l2004l'l. 2005hk 

(IChornock et al.j 20061: Sahu et al.l l2008l: 

McCullv et al.l 

2014all. 2008A ImcCuIIv et al. 2014al'l. 2008 

la and 2010ae 
gll. Sill, SII, 

iStritzinffer et al.l l2014h: Ol, Oil, Nal, M 


Call, Till, Fell, Felll and Coll. Not all ions were detected 
in all spectra, e.g. the ones with higher ionization poten¬ 
tial (On and Felll) were not found in the late-phase -|-17d 
spectrum. 

For both models, the optimum set of parameters that 
give the best fit for the synthesized spectrum to the ob¬ 
served one were found automatically by SYNAPPS via x^~ 
minimization. The total number of the optimized parame¬ 
ters were ~ 50 in a typical run. Note that there is an es- 
sential degeneracy betw e en T(vref), v-o hot (or Vmin) and Ve 
d Jeffery fc BrancM 1 19901 : [P arren ^ l2014i ~). thus, their finally 
adopted values may not represent a unique solution for any 
spectra. 

The final best-fit model spectra for both Model A and 
Model B are plotted together with the observations in Fig.[Sl 
Fig. [S]shows the single-ion contributions to the best-fit mod¬ 
els for the -bOd spectrum, i.e. at the phase of V-band maxi¬ 
mum. The basic global parameters for both models are col¬ 
lected in Table [S] 

From Tableland Figs [5]-[6] it is immediately apparent 
that, despite having quite different Vphot values, both Model 
A and B give acceptable hts to the observed spectra. Both 
models reproduce the major spectral features well, although 
neither fits perfectly all the weak, narrow humps that be¬ 
come stronger in the post-maximum spectra. This is some¬ 
what surprising, because the two models converged to Vphot 
values that are different by ~ 3000 km s“^. Since the uncer¬ 
tainty of the photospheric velocities derived by SYNAPPS 
are thought to be ~ 500 - 1000 km s“^, at least for SNe la 
llParrentjboiJl ). such a high level of ambiguity, at first, seems 
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Figure 5. The best-fit SYNAPPS model spectra for Model A (solid line) and Model B (dashed line) plotted together with the observed 
spectra at —9d, —2d, Od, -|-3d and -|-17d phases. Strong blending makes feature identification ambiguous at all phases. Some of the 
strongest contributing ions are marked by vertical lines and labels. 

to be unexpected. It is, how ever, not unprecedented among 
other SNe lax: for example, IStritzinger et al] (l2014l) found 
similar inconsistency between their SYNAPPS models for 
the simulteneous optical {vphot ~ 4200 km s '^) and near-IR 
{vphot ~ 1700 km s”*^) spectra of SN 2010ae at -|-18d phase. 

A possible (and likely) explanation for this issue is that 
all modeled SN 2011ay spectra consist of heavily blended 
features, even at very early phases. There are no individ¬ 
ual, unblended, single-ion features visible that can be iden¬ 
tified unambiguously, except maybe OI A7775 in the pre¬ 
maximum spectra and the Call near-IR triplet after maxi¬ 
mum. In both models considered here the presence of iron 
(both Fell and Felll) can be found at —9d, and blend¬ 
ing with Fell and Coll becomes excessively dominant at 
after maximum light. As noted bv iBranch et al.l 
“iron curtain” of strong Fell prevented the 
secure identification of many other features in the post¬ 
maximum spectra of SN 2002cx, and this is exactly what we 
found here for SN 2011ay. The projected Doppler-velocities 


later p hases, 
(120041) ■ the 


of blended features become quickly ill-constrained with the 
increasing number of overlapping lines from different ions. 

The evolution of the minimum velocities of the line 
forming region for each species in Model B are plotted in 
Fig. [3 On the contrary. Model A assumed that all features 
are formed with a minimum velocity of Vphot as given in 
Table [S] It is apparent that while in Model B the photo- 
spheric velocity is significantly lower than that in Model A, 
Model B contains some low-mass ions (OII, Nal, Mg II) that 
are formed at higher velocities. The minimum velocities of 
these “detached” features are close to Vphot of Model A, but 
they tend to decrease below that, and converge toward the 
photosphere at later phase s (except for Oil, w hich stays 
constantly at high velocity). iBranch et akl ||2004|) also found 
similar “high-velocity” features in the spectra of SN 2002cx, 
but only at post-maximum phases. Fig. [7]suggests that such 
detached features forming at higher velocities may have ap¬ 
peared in SN 2011ay as early as 9-10 days before maximum. 

From Fig. [5] and [6] it is seen that low-Z elements, such as 
Ol, Nal, MgII, start appearing at very early phases (more 
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Rest wavelength (A) 


Figure 6. Single-ion contributions to the Od spectrum fit (solid line = Model A, dashed line = Model B). 


than one week pre-maximum), and they remain detectable 
even at post-maximum. SII may also be present (at least in 
Model B), which is often cited as a sign of thermonuclear 
explosion. These features are frequently observed in SNe la, 
but they are absent in SNe Ib/c. 

iFolev et al.l (l2013h reported the possible presence of car¬ 
bon in SN 201 lay via the detection of CII A6580 and A7234 
features, e ven though these f eatures appear very weak (see 
Fig. 23 in iFolev et al.l [2OI3I I. We have checked the pres¬ 
ence/absence of C I, CII or C III using SYNAPPS. The —9d 
and —2d spectra were selected for this purpose, because car¬ 
bon features are expected to be the strongest during the 
pre-maximum phases. For this test we adopted Model B, 
be cause that model has similar velocities to those adopted 
by IFolev et al.l (l2013fl . 

As a first step, we simply added either Cl, CII or CIII 
to the best-fit Model B spectra to get carbon-enhanced mod¬ 
els. The purpose of this step was to check which carbon 
features should be visible under the physical conditions cor¬ 


responding to Model B. Fig.[8]shows these carbon-enhanced 
models together with the carbon-free models and the —2d 
and —9d observed spectra. Note that the carbon features are 
artificially enhanced for making their identification easier. 

It is seen that large amount of CI or C III would cause 
observable features only at wavelengths longer than 8000 A 
where the observed spectra have the lowest signal-to-noise 
and contamination from telluric lines is the strongest. Thus, 
the reliable identification of these ions is not possible from 
these spectra. 

However, the CII A6580 and A7234 features offer better 
opportunity to detect carbon in the optical spectra. There is 
a feature in both the —9d and the —2d observed spect ra that 
could be CII A7234, as found bv IFolev et al.l ll2013li . Note 
that this feature is originally modeled with high-velocity OII 
in both Model A and B, which is not a secure identification 
(see below). Even though the simultaneous presence of the 
CII A6580 feature is not observed in the —2d spectrum, there 
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Rest wavelength (A) 


Figure 7. Evolution of the minimum velocities of the line-forming 
regions for each species in Model B. The uncertainty of these ve¬ 
locities is at least ~ 1,000 km s“^. The lowest velocity at each 
phase corresponds to the photospheric velocity at the given epoch. 
For comparison, the photospheric velocity in Model A is also plot¬ 
ted with asterisks. 


Figure 8. Modified Model B spectra containing artificially en¬ 
hanced carbon features (thick solid and dashed lines) together 
with the —2d observed spectrum (grey curve) and the original, 
best-fit Model B spectrum without carbon (thin solid line). The 
observed and carbon-free model spectra at —9d are also shown 
for comparison. 


may be a weak feature at the expected position in the —9d 
spectrum. 

In order to test whether it is more appropriate to ex¬ 
plain the feature around 7100 A with CII rather than Oil, 
we replaced OII with CII in the best-fit — 9d Model B 
spectrum, and recomputed the fitting. The excitation tem¬ 
perature for carbon was set variable to take into account 
the possibility for non-LTE excitations. It was found that 
Tea;c = 20, 000 K (which is a factor of 2 higher than Tphot at 
this phase) would indeed reduce the strength of the A6580 
feature relative to that of the A7234 one. 

In Fig. [9] the new carbon-enhanced model is compared 
with the original best-fit Model B spectrum (without car¬ 
bon) and the observations. It is seen that while the CII 
A7234 feature is a similarly good fit to the observed spec¬ 
trum as the high-velocity OII feature in the original model, 
this is not true for the CII A6580 feature, the latter be¬ 
ing much stronger than the observed notch around 6400 A. 
Note that setting the excitation temperature of CII close 
to the photospheric temperature, Tphot = 10, 000 K, would 
make the A6580 feature even stronger, enhancing the incon¬ 
sistency between the carbon-enhanced model and the ob¬ 
servations. Thus, if CII were responsible for the observed 
feature at ~ 7100 A then the A6580 feature should appear 
much stronger than observed. 

It is concluded that our SYNAPPS modeling does not 
confirm the detection of CII features in the pre-maximum 
optical spectra of SN 2011ay, and the same is true for both 
Cl and cm. 

The identification of the A7100 feature as high-velocity 
On in the early-phase spectra (see Fig. [^and[7]l is also un¬ 
certain. OII is usually absent in the pre-maximum spectra 
of SNe la, but it was found in other SNe la x around maxi- 
mum light, namely in 2008A and 2005hk, bv lMcCullv et akl 
ll2014ah . Given that our models for SN 2011ay have photo¬ 
spheric temperatures not exceeding 10,000 K (cf. Table [5]), 
the appearance of O II with nearly the same str ength as OI is 
not ex pected, at least in quasi-LTE conditions. iHatano et al.l 
lll999l) predicted equal optical depth for these two ions at 



Rest wavelength (A) 


Figure 9. Comparison of the modified Model B containing C II 
instead of high-velocity Oil (dashed line) with the original Model 
B (having Oil, without CII, thick solid line) and the observed —9d 
spectrum (grey curve). If CII A7234 were present, A6580 would 
appear much stronger than observed (see text for details). 

T ^ 12, 000 K, but note that they used tie ~ 5 x 10® cm“^ 
for the electron density, a typical value for SNe la-s. If the 
electron density in SN 2011ay were lower than this, then ion¬ 
ization would be more effective at a given temperature, thus, 
the OII features could become stronger even at T < 10, 000 
K. Nevertheless, since it is based on only a single feature, 
the presence of Oil in SN 2011ay cannot be proven unam¬ 
biguously. 

Regarding the intermediate-mass elements (IME), Call 
is close to the detection limit in the pre-maximum spectra, 
but its near-IR triplet (IR3) A8579 feature becomes very 
strong by two weeks after maximum. Similar behavior is ob¬ 
served in overluminuous S Ne Ia-91T (see e.g. IPhillips et akl 
I1992I : [Caravini et al.ll20()3 ). but those objects are character¬ 
ized by very different light curve properties. 

This kind of spectral evolution is the opposite of wha t 
is observed in normal SNe la (see e.g. iMarion et akll^Ma ): 
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in most SNe la the Ca II IRS starts as a strong high-velocity 
feature (HVF) at n ^ 20,000 km s“^ at about two weeks 
before maximum, while low-Z elements are absent or very 
weak. By the time of maximum light the Call HVF disap¬ 
pears, while its photospheric component (PVF, v ~ 10, 000 
km s~ ) gets stronger. OI A7775 ca n, sometimes, be de - 
tected as a HVF at very early phases llParrent et al.l[2012ll . 
but it is frequently absent and its PVF shows up only after 
maximum light. 

It is seen in the spectra presented in this paper that 
SN 2011ay does not exhibit such kind of HVF for any ion 
at any phase observed. This also appears true for all other 
SNe lax observed so far llFolev et al.ll201§ l. It is interesting 
that the absence of Call HVF is also characteristic for the 
low-velocity, 91bg sub type of SNe la llChildress et al.|[20l4 
ISilverman et al.ll2015t l. 

Unlike the Call features, the Fell lines appear quite 
strong in SN 2011ay, even during pre-maximum. The strong 
presence of Fell at such early phases were not observed in 
most “normal” SNe la spectra. Also, they were not iden- 
tifled in other pre- maximum lax spectra, e.g. SN 2002cx 
llBranch et ahlliooll or SN 2005hk llSahu et al.il200^ 1. Their 
strong appearance in SN 2011ay might be due to its cooler 
photospheric temperature during pre-maximum: Tphot ^ 
10,000 K (Table [5]). SNe la typically have Tphot 10,000 
K before maximum. Concerning other SNe lax. Fell fea¬ 
tures were fo und strong in the post-maximum spectra of 
SNe 2002cx (iBranch et al.l l2004l l and 2005hk llSahu et al.l 
l2008ll . when the photospheric temperature cooled below 
~ 10,000 K. After maximum light, Fell and Call features 
dominate the observed spectra of SN 2011ay (Fig. [5]), similar 
to SNe la. 

3.3 Measuring velocities in SNe lax spectra 

One of the most intriguing properties of SNe lax is their 
low e xpansion velocities fe.g. lBranch et aLll2004l : rFolev et al.l 
l2013h . which can, in extreme case s, be as low as ~ 
2000 km s“^ (IStritzinger et all 120141 '). This has been con¬ 
firmed for multiple SNe lax by spectral modeling with 
SYNOW/SYNAPPS, i.e. similar methodology to that which 
has been applied in this paper. To esti mate the expansion 
velocity for a bigger sample of SNe lax, iFolev et al.l (l2013l l 
used the projected Doppler-shift of the absorption minimum 
of the feature around A6200 A, which was interpreted as 
due to Sill A6355, similar to Type la SNe. Based on this 
assumption, iFolev et al.l ll2013ll concluded that the “ejecta 
velocity”, |u|, defined this way is less than 8000 km s“^ for 
all SNe lax. All SNe lax spectra studied to date indeed show 
signs for significantly lower expansion velocities than most 
SNe la which usual ly have usiii > 10, 0 00 km s“^ around 
maximum light (e.g. iBlondin et ^l2006h . however, we show 
below that such “quick-look” velocity estimates based on 
single spectral features may lead to significantly under- or 
overestimated velocities for SNe lax, which have much more 
complex spectra than SNe la. 

We focus on the feature appearing around ~ A6200 A 
in the pre-maximum spectra of SN 2011ay. Fig. HOI shows 
this feature before and around maximum light compared to 
two other SNe lax observed at similar phases. It is seen 
that this feature is broad and shallow. In this respect, 
SN 2011ay and some other SNe lax look similar to SNe la- 


Si 11X6355 velocity (10^ km s“') 

-16 -10 -5 0 5 



Rest wavelength (A) 


Figure 10. Comparison of the near-maximum spectra of three 
different SNe lax in the region of 6000-6500 A: SN 2005hk 
llPhillips et al.l 1200 7f) . SN 2 008A t|Blon din_gt_alJ |2012|1 . and 
SN 2011ay (present work and IFolev et al. I I2013I) . All epochs are 
given relative to V maxima. The spe ctra of SNe 20Q5hk and 2Q 08A 
were downloaded from WISeREP dYaron &: Gal-YamI l20l3l . A 
velocity scale relative to Si II A6355 is also shown on the upper 
horizontal axis. 


91T llBranch et al.ll20o3 : IFolev et al.ll2013l l. i.e. they usually 
show weaker Sill A6355 lines compared to normal la-s. How¬ 
ever, the post-maximum spectral evolution of the two classes 
of SNe are markedly different. 

The first issue with the velocity measurement from 
this fe ature is related to its line strength. Ijefferv fc Brancbl 
lll99dl studied the effects of line strength, ejecta density and 
other parameters on the Doppler-shift of the absorption min¬ 
imum of P Cygni features in SN atmospheres. They showed 
that for features having optical depth r ~ 1 at the photo¬ 
sphere, the absorption minimum forms very close to Vphot- 
That is why the velocity from the Si II A6355 feature gives so 
reasonable photospheric velocities in SNe la: this feature is 
stron g and relatively unblended in SNe la (e.g. lParrent et ahl 
I2OI2II . For weak features having r < 1, however, resonant 
scattering tends to dominate over absorption, which may 
shift the absorption minimum significantly below Vphot- On 
the contrary, much stronger features, like Ca II, tend to show 
absorption minima at velocities much higher than Vphot, as 
expected. 

Fig. [To] suggests that SN 2011ay had similar photo¬ 
spheric velocity to SN 2008A, and both of them expanded 
faster than SN 2005hk. Using SYNAPPS, iMcCullv et Zl 
ll2014ah derived ~ 7000 and ^ 8500 krn s~^ for 2005hk 
and 2008A, respectively, while l^lev et all (l2013ll estimated 
~ 4500 and ~ 6400 km s“^ from the A6200 feature. The 
~ 2000 km s“^ difference between the single-feature and the 
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Si II W365 velocity (10^ km s“') 

-16 -10 -5 0 5 



Rest wavelength (A) 


Figure 11. Spectral evolution of SN 201 lay between 6000-6500 
A. The observed spectrum (thick gray line) is plotted together 
with the SYNAPPS model from Model A (solid line) as well as 
single-ion contributions from Si II (dashed line) and Fe II (dot¬ 
ted line). Fe II A6456 starts blending with Si II A6355 a few 
days before maximum, and it becomes the dominant feature post¬ 
maximum. A velocity scale relative to Si II A6355 is also shown 
on the upper horizontal axis. 


spectrum-modeling velocity illustrates the issue of the veloc¬ 
ity measurement in SNe lax spectra. For SN 2011ay we find 
a similar conflict: at maximum light the single-feature ve¬ 
locity estimate g ives ~ 5600 km s“^ (ISilverman et al.ll201ll : 
iFolev et'aDl2013 ). while our Model A predicts Vphot 9000 
km s“^ (with an uncertainty of at least 1000 km s“^). 
It seems that SYNAPPS modeling tend to produce higher 
velocities than the single-feature “quick-look” estimates. 

However, there is an even bigger issue with the velocity 
measurement in SNe lax spectra. This is the severe blending 
across the whole optical regime, which prevents any unique, 
secure line identification, as was detailed in the previous 
section. This affects the A6200 feature, further plaguing the 
velocity determination. 

Fig. [TT] again zooms in the 6000 - 6500 A regime while 


Si II k6365 velocity (10^ km s“') 

-15 -10 -5 0 5 



Rest wavelength (A) 


Figure 12. The same as Fig. llll but showing Model B, and adding 
the single-ion contributions from S II and Co II. Complex blend¬ 
ing between these features results in nearly the same observed 
spectra, even though the photospheric velocities are ~ 3000 km 
s“^ less than in Model A. 


showing the observed spectral evolution of SN 2011ay to¬ 
gether with the model spectra by Model A (i.e. the one 
having higher Vphot)- From the single-ion contributions it 
is seen that the broad A6200 feature is due to blend¬ 
ing between Sill A6355 and Fe II A6456, even during 
the pre-maximum phases. Similar blending is also visible 


in other SNe lax ^ as found in 2002cx at -|-7d {Tbb ~ 
9000 K, iLi et al.l l2003l : Branch et al.l 20041). 2005h k at - 
Id {Tbb ~ 90 00 K, [milips et al.ll2007l : ISahu et'Zll2008l) . 
2008A at -3d jM cCully et al.l 2014a ). and 2010ae at -|-16d 
llStritzinger et al.l l2014l ). This blending broadens the ob¬ 
served profile and shifts the middle of the feature toward 
lower velocities, which may likely explain the discrepancies 
between the “quick-look” and the modeled velocities men¬ 
tioned above. 

Moreover, this is not yet the full story, because, as il¬ 
lustrated in Fig. 1121 blending between the different single¬ 
ion features in this regime may result in nearly the same 
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observed spectrum, even when the photospheric velocity is 
close to ~ 6000 km s“^, i.e. close to what the “quick-look” 
velocity estimate predicts. Thus, even when using a sophis¬ 
ticated parametrized spectrum modeling code and putting 
in all possible features that likely contribute to the observed 
spectrum, the resulting model may still be ambiguous: if the 
observed spectrum consists of broad features and almost all 
of them are complex blends, then the photospheric velocity 
becomes rather ill-constrained. 

A good example for this ambiguity is SN 2011ay: as 
shown in the previous section, its observed spectra can be 
explained either with Vphot 6000 or ~ 9000 km s“^ al¬ 
most equally well. This and the appearance of the detached 
features (those having higher velocities than Vphot) in the 
lower velocity Model B might also be a hint for a multiple 
velocity structure in the ejecta, i.e. a disk-like or a jet-like 
configuration. Investigations of such models are beyond the 
scope of this paper, and might be possible only by having 
more extensive data coverage, both in wavelength (I.e. by 
adding infrared and/or UV spectra) and in time. 


3.4 Spectral energy distributions and light-curve 
modeling 

To calculate the SEDs of SN 2011ay, we converted the BVRI 
and Swift /UVOT magnitudes to flux es using the cali¬ 
bration of iBessell. Castelli fc Pl^ (Il998l l. The fluxes were 
der eddened using the ga l actic reddening law parametrized 
by [Fitzpatrick &: Massal ll200lf l assuming Rv = 3.1 and 
adopting E{B — V) — 0.081 mag (§3.1.). The combined UV- 
optical SEDs are compared with the three HET spectra in 
Figure [T3l After correcting for the small uncertainties in the 
absolute flux calibration of HET data, the spectra and the 
broadband SEDs match very well. Trustworthiness of Swift 
data at -1-17 days are lower, because SN 2011ay was already 
quite faint in the UV at that time. 

The quasi-bolometric light curve was derived by in¬ 
tegrating the dereddened values of the combined UV- 
optical SEDs against wavelength. The long-wavelength 
contribution (not covered by these data) was estimated 
by htting a Rayleigh-Jeans tail to the red end of the 
observed SEDs, and integrating it to infinity. In ad¬ 
dition, the light curve was supplemented by two pub¬ 
lished KAIT-measurements: a non-detection on March 1 
(2 455 6 21.7 JD), and the unfi ltered magnitude at the dis¬ 
covery telanchard et al]|201l[l . The luminosity calculated 
from the latter was used as a lower limit at that epoch. Fi¬ 
nally, the Integrated fluxes were converted to luminosities 
using D = 86.9 Mpc (see §1). 

We appl ied a generalized analytic light curve mode l 
published by Chatzopoulos. Wheeler fc Vinkdl ll2009l . l2012l l 
to derive the main parameters of the SN. This model i s based 
on the radioactive decay dif fusion model of ArnettI (Il982l l 
and its generalized form by IValenti et al.l (l2008l l. but also 
takes into account the gamma-ray leakage from the ejecta. 
In this model the output luminosity can be expressed as 


L{t) = Mnic 


(cNi — eco) / 22;e^ dz + eco 


X / 2ze' 


2"^ —2yz-\-2za 


dz 


(3) 


Table 6. Physical parameters determined from the modeling of 
the bolometric light curve. 


to (JD) 

2 455 633.0 ± 1.5 

Date of explosion 

^rise (d) 

14 ± 1 

Rise time to maximum 

MnI (Mq) 

0.225 ± 0.010 

Initial ^®Ni mass 

Aa 

-^norm 

25 ± 7 

7 -ray leakage parameter 


Notes. “The parameter connected to gamma-ray leaking (see Eq. 
[3|and text for details). 


where x = t/tm, tm is the effective diffusion time (roughly 
equal to the rise time to maximum), y = tm/{2tm) with 
tNi = 8.8 days, s = tm(fco - ^Ni)/(2tcoiNi) with tco = 111.3 
days, Mni is the initial nickel mass, eni = 3.9 x 10^'^ 
ergs^^g”’^ and eco = 6.8 x 10® ergs“^g“^ are the energy 
generation rates due to Ni- and Co-decay. 

The last term of Eq. describes the amount of 
gamma-ray leaking. Assuming a spherical uniform density 
ejecta with radius R = vt, homologous expansion (due to 

scaling), and the Ni/Co confined in the center, yields 
A = (3«:-yMej) / where is the gamma-ray opacity, 

Mej is the ejecta mass. Note that we used the dimensionless 
form of (l — exp (—Anorm • tfo)), where Ho = t/10 days. 

The quasi-bolometric light curve and the best-htting 
model are shown in Figure 1141 The corresponding parame¬ 
ters are presented in Table |6l Both the nickel mass (~ 0.22 
Mq) and the rise time (~ 14 d) are lower than the usual 
values for SNe la (~ 0.6 M q and ~ 18 d, re spectively) 


similar to the conclusion by Folev et al.l (l2013l l. The short 
rise tim e implies low eject a mass. Following the method ap¬ 
plied bv iFolev et al.l (120091 1 for SN 2008ha and assuming that 
the ejecta of SN 2011ay has the same mean opacity as that 
of a normal SN la, we have 


V Ula / y^rise,Ia J 


(4) 


and 


Ula \ Crise.Ia 


(5) 


Adopting V ~ 9,300 km s“^, as suggested by the higher- 
velocity Model A (§3.2) and trise = 14 days together with 
the reference values of Ekin.ia = 10^^ erg, Mej,ia = 1.4 Mq, 
Via, = 10 000 km s”*^ dFolev et al.ll200 9|l and the average value 
of tri se.ia = 18.0 days (see Ganeshalingam. Li fc Filippenkol 
these formulae result in Ekin = 0.5-10®'^ erg and Mej 
= 0.8 M© as the kinetic energy and the ejecta mass of 
SN 2011ay, respectively. Alternatively, using v ~ 6000 km 
s“^ from the lower-velocity Model B results in Mej ~ 0.5M© 
and Ekin = O.ITO®^ erg. The ejecta mass and expansion en¬ 
ergy estimated this way are consistent with previous find¬ 
ings for SNe lax. The lower nickel mass, Mni ~ 0.2 M© as 
calculated from the Arnett-model, also suggests that the ex¬ 
plosion, if indeed thermonuclear, may be less energetic that 
in Type la-s. 
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Figure 13. Scaled HET spectra and the UV-optical SEDs of SN 2011ay at different epochs relative to V maximum: -9 days (2011 March 
24, top left), -4 days (2011 March 29, bottom left), +17 days (2011 April 19, bottom right). 



JD - 2 450 000 

Figure 14. Quasi-bolometric light curve of SN 2011ay with the 
best-fitting Arnett-model (solid line) and the one with assuming 
full gamma-ray trapping (dashed line). The fir st two luminosi¬ 
ties w ere calculated from the measurements of iBlanchard et al.l 
(l201lh . see the text for details. 


4 DISCUSSION AND CONCLUSIONS 

Based on photometric and spectroscopic data obtained dur¬ 
ing the early phases we carried out a detailed analysis on 
SN 2011ay, one of the members of the recently defined class 
of SNe lax. The spectra as well as the light and color curves 
are similar to those of other objects belonging to this group 
(see lFolev et al.l[2013l . and references therein). 

We calculated model spectra with the parametrized 


modeling code SYNAPPS to get more information on the 
spectral evolution and the physical properties of the ejecta. 
As presented in Section 13.2[ the spectral characteristics of 
SN 2011ay are basically similar to those of other SNe lax. 
We found that Fell features appear even before maximum, 
which can be explained with the relatively lower photo- 
spheric temperature {Tphot ~ 8,000 K) of the ejecta. 

The presence of Fell and other features that cause se¬ 
vere blending across the entire visible spectral domain makes 
the direct analysis of such spectra very difficult and ambigu¬ 
ous. We found that it is possible to fit most of the broad 
features of SN 2011ay by two models (Model A and B) that 
have similar chemical composition but very different photo- 
spheric velocity (~ 9,300 km s”*^ and ~ 6,000 km s“^ at 
maximum, respectively). 

The effect of strong blending also makes the “quick- 
look” velocity estimates uncertain. We showed in Section 
m\ that measuring the absorption minimum of the A6200 
feature and assuming that it is due to purely Sill A6355 
might underestimate VpUot, because of the effect of blending 
with Fell and Coll. This does not question the fact that 
SNe lax generally have lower expansion velocities than SNe 
la, but the non-uniqueness of the spectrum modeling fits 
presented in this paper makes the velocities of such events 
rather ill-constrained. 

Fitting a radiative-diffusion model (taking into account 
gamma-ray leaking) to the quasi-bolometric light curve of 
SN 2011ay resulted in physical parameters for the rise time 
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to maximum, ejecta mass, initial nickel mass and kinetic 
energy that are all below their me an value for SNe la, as 
pointed out by I Foie v et alJ (l201.ll ~). Although SN 2011ay 
seems to be close to the upper limits of these parameters, 
it clearly belongs to the recently defined class: the values 
of ®®Ni mass and total ejecta mass (0.225±0.010 M© and 
Mq, respectively) are similar to the parameters of 
other, brigh t (My > -18 mag) SNe lax (SNe 20 05hk, 2008A, 
2009ku; see I Foie v etliD l2013l : iMcCullv et alJ 12014^ . Not e 
that in addition to the criteria defined by Folev et alJ ll2013ll . 
the faster decline rate of the bolometric light curve, indi¬ 
cated by the relatively low value of the gamma-ray trapping 
parameter Anorm (see Table [6|, may also be characteristic 
for SNe lax. 

The results of the light curve modeling strengthen the 
picture that the progenitor of SN 2011ay is likely an incom¬ 
pletely exploded white dwarf, which seems to be the best 
explanation of this kind of SNe to date. The presence of 
Sill and SII features in the early spectra of SN 2011ay in 
Figure [5] may be signs of thermonuclear explosion, which 
might give additional support for the white dwarf explosion 
scenario. 

These results, as well as other detailed studies of single 
events, may help us to understand better the properties of 
type lax SNe. However, there are still a lot of open ques¬ 
tions about the true nature and origin of this class of stellar 
explosions. 
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